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Introduction 


Out of the welter of orbital data from a satellite it is frequently 
useful to be able to isolate and examine the behavior of one phenomena. 

We have been particularly interested in knowing where perigee is located 
at any particular time in the lifetime of the 0G0-2 satellite, particularly 
as a function of local time. 

There are several reasons for this. First, the orbital plane of a 
polar satellite moves very slowly with respect to the earth sun line, i.e. 
the satellite position varies very slowly in local time. This fact is 
significant when attempting to determine diurnal effects (as to eliminate 
such effects) from experimental satellite data. Similarly, most geophysical 
phenomena are latitude dependent and some phenomena may only be observed 
at low or high latitudes. (For example, the equatorial electro jet phenomena 
in magnetic field studies is both local time and latitude dependent). It 
is therefore important for the experimenter to take orbital factors into 
account in planning his experiment and in analyzing his results. 

If a satellite were in orbit about a perfectly spherical earth with 
no outside perturbations, its perigee would maintain a fixed latitude and 
the local time of the orbital plane would depend only on the right ascension 
of the sun. In fact, however, the earth is not spherical and the orbit is 
subject to many outside perturbations. 

As a result of these perturbations, the orbital elements are no longer 
constant. By considering a first approximation of those orbital changes 
which have the major effects on the latitude and local time of perigee, we 
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have devised a simple calculation giving results which, for a nine-month 
period, check to within one degree in both latitude and local time with 
a more sophisticated computation. 

Method 

If we consider three coordinate systems (Figure 1) such that X, Y, Z 
are the geocentric inertial system of some epoch with X along the equinox, 
X 1 , Y 1 , Z 1 are in the plane of the orbit with X 1 toward perigee and Z 1 
chosen so that its angle with Z is less than 90°, and X 2 , Y 3 , Z 2 are a 
rotation of X, Y, Z so that X 3 lies along the projection of the earth-sun 
line in the X, Y plane. 

Where : 

Q = right ascension of ascending node of the satellite 
ao = argument of perigee 
i = inclination 

RA = right ascension of the mean sun 

The transformation from the X 1 system to the X system consists of 
three rotations 1) through uo, 2) through i, and 3) through Q. The trans- 
formation from the X to the X 3 system is a rotation through RA. The 
resulting transformation matrix from the X 1 to the X 2 system is then given 
by: 
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Where: 


A 11 
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cos (RA 

- n) + 

sin(tu) 

cos (i) 

sin(RA 
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A 12 
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cos (RA 
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sin(RA 

- Cl) 

A 13 

= 

- sin(i) 
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- sin(i) 
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sin(i) 

sin(u)) 
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Now if a is the earth central distance to the satellite at perigee 
then in the Y 1 system perigee is at 


Therefore, applying the matrix A, perigee in the X 2 system is 


( 2 ) 


a [Z COS ^ C0S ^ A - Q) - sin (jd cos(i) sin(RA - Q)^| 
a[£cos(aj) sin(RA - Q) + sin(uj) cos(i) cos(RA - Q) j 
a Q sin(i) sin(u)7] 

Now in the X 2 system the tangent of (local time - 180°) equals the X 2 
coordinate of perigee divided by the Z 2 coordinate of perigee. 


or 

TAN (L.T. 


180°) = sin (uu.) cos i cos(RA - Q) - cos (D sin(RA - Q) 
sin(uu) cos i sin(RA - Q) + cos (D cos(RA - fi) 


Also SIN (Latitude) = 
and 


coordinate of perigee 


sin i sin a) 


(3) 

(4) 
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TAN(Latitude) 
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With the above results, then, we may easily and exactly compute both 
the latitude and local time of a satellite if we know the values of Q, ub 
i, RA, and a. Normally, however, one is given these values at some epoch, 
along with values of u) and fi. We have then used the following approxima- 

9 

tions. If we are given u) 0 , ftp, fig, u)g , i, RA 0 , and a at some initial 
time, t 0 , we have assumed i, a, iOg , and C^ to remain constant and have 
computed Q, tu, and RA at time t using: 

RA = 360/365.256 

Cl = Qq d" Clg (t “ tg) 

( 6 ) 

CU = (JJg + U)g (t - tg) 

RA = RAg + RAg (t - tg) 
and then substituting in (3) and (5) above. 


• » 

If Cz and u) are unknown, they may be approximated readily. 
From (1) page 32: 


Where: 
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a s (1 - e 2 )^ 
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.606546 
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Universal gravitation constant 
Mass of the earth 
Semi -major axis of orbit 
Eccentricity 
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Inserting the constants gives 

3kN = 20.8158 (10 13 ) 

a 2 (1 - e 2 ) 2 a 3,5 (l - e 2 ) 2 

where a is in kilometers. 

Discussion 

Conrath (Reference 1) classifies perturbations as either secular or 
periodic and states (Pg. 12): "In considering the behavior of the orbit 

as a whole, the secular perturbations provides the best picture, and the 
periodic terms can be neglected. However, if the calculations of an 
ephemeris is desired, the periodic perturbations must be considered." Since 
we wish to construct a picture of how the orbit is situated with respect to 
local time, we have considered that a first approximation to the secular 
terms (a) and Q) is sufficiently accurate. These two variations primarily 
result from the R" 3 term in the potential expansion of the earth's gravita- 
tional field and, therefore, from the non-sphericity of the earth. From 
the approximate expressions given in (7) it may be seen that the position 
of the ascending node regresses in proportion to the cosine of the inclina- 
tions, and is therefore quite small for a polar satellite, and that the 
motion of arguments of perigee depends upon the sign of (2-2.5 sin 2 i). It 
is thus constant for i = 63° 26' . Plots of the rate of change of these 
angles of a function of i may be found in Reference 1, pages 42-48. 

For our computations we have used the right ascension of the (fictitious) 
mean sun (Ref. 4, Pg. 139-141). This brings our computations into agreement 
with the computation of Universal Time where the mean sun is also used 
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(Ref. 3, Pg. 139-141, Ref. 4, Pg. 73-76). The actual sun may differ from 
the mean sun by a maximum of approximately 3° or 12 minutes of local time. 

In order to compute the "actual" local time the apparent right ascension 
of the (real) sun should be substituted for our right ascension of the 
mean sun. The apparent right ascension of the real sun may be found in 
the American Ephemeris and Nautical Almanac, Pages 18-33; the right ascen- 
sion of the mean sun is equivalent to the longitude of the mean sun (Ref, 3, 
Page 140) tabulated on Pg. 50 of the Almanac. 

Results 

The 0G0-2 satellite was launched October 14, 1965 into a low polar 
orbit with the following characteristics (0G0-2 operations center, Nov. 23, 
1965) at epoch October 24th, 0 hrs , , 0 min., 0 sec. : 

a = 7340.5 km 

i = 87.359 degrees 

RA 0 = 208.26 degrees 

Qo = 280.49 degrees 

U0o = 144.211 degrees 

Q = - . 2839°/ day 

(L = -3 . 0476°/ day 

Using these figures, we have produced the plot shown in Figure 2. 

Each dot represents perigee position at zero hours on a particular day 
with the first day of each month labeled. The black dots indicate perigee 
in the southern hemisphere and the open dots indicate perigee in the 
northern hemisphere. The plot is polar with 9 representing local time 
and r representing latitude. A comparison of our results for the 0G0-2 
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satellite with the ORB-1 printout from the GSFC Theory and Analysis Office 
showed agreement in both latitude and local time within one degree for a 
period of nine months. 

Summary and Conclusions 

We have developed a procedure for a graphic display of one specialized 
feature of a satellite orbit, the relationship of the perigee of the orbit 
to latitude and to local time. From such a plot (e.g. Fig. 2) it is a very 
simple matter to see exactly when in a satellite lifetime it will encounter 
a particular sunlight configuration. One may also readily see at what 
latitude the satellite approaches closest to the earth. Obvious applications 
are predictions of eclipse times and periods of maximum drag. From a geo- 
physical viewpoint, one may predict when satellite data should be examined 
for effects which are local time and/or latitude dependent. 
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